Zero-bias anomaly in a nanowire quantum dot coupled to superconductors 
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We studied the low-energy states of spin- 1/2 quantum dots defined in InAs/InP nanowires and 
coupled to aluminium superconducting leads. By varying the superconducting gap, A, with a 
magnetic field, B, we investigated the transition from strong coupling, A <C Tk, to weak coupling, 
A ^ Tk, where Tk is the Kondo temperature. Below the critical field, we observe a persisting 
zero-bias Kondo resonance that vanishes only for low B or higher temperatures, leaving the room to 
more robust sub-gap structures at bias voltages between A and 2 A. For strong and approximately 
symmetric tunnel couplings, a Josephson supercurrent is observed in addition to the Kondo peak. 
We ascribe the coexistence of a Kondo resonance and a superconducting gap to a significant density 
of intra-gap quasiparticle states, and the finite-bias sub-gap structures to tunneling through Shiba 
states. Our results, supported by numerical calculations, own relevance also in relation to tunnel- 
spectroscopy experiments aiming at the observation of Majorana fermions in hybrid nanostructures. 



Hybrid devices which couple superconducting (S) elec- 
trical leads to low-dimensional semiconductors have re- 
ceived great attention due to their fascinating underlying 
physics Further interest in this field has been gener- 
ated by recent theoretical predictions on the existence of 
Majorana fermions at the edges of one-dimensional semi- 
conductor nanowires (NWs) with strong spin-orbit inter- 
action connected to S electrodes [1]. Zero-bias conduc- 
tance peaks meeting some of the expected characteristic 
signatures of Majorana physics were recently reported in 
hybrid devices based on InSb [1, 01 and, more recently, 
InAs NWs [Hi- In the past years, quantum dots (QDs) 
coupled to superconducting leads have been widely ex- 
plored as tunable Josephson junctionsP, 0, or as build- 
ing blocks of Cooper-pair splitters i-e. devices that 
could serve as sources of entangled pairs of electrons. Hy- 
brid superconductor-QD devices also constitute versatile 
platforms for studying fundamental issues, such as the 
physics of the Andreev bound states (ABS) [13, [HI or 
the interplay between the Kondo effect and the super- 
conducting proximity effect |l^ - l2]| . 

The Kondo effect usually stems from the antiferromag- 
netic coupling of a localized electron spin and a Fermi sea 
of conduction electrons. Below a characteristic tempera- 
ture Tk, the so-called Kondo temperature, a many-body 
spin-singlet state is formed, leading to the partial or com- 
plete screening of the local magnetic moment. This phe- 
nomenon, discovered in metals containing diluted mag- 
netic impurities, is now routinely found in individual QDs 
with a spin-degenerate ground state, e.g. QDs hosting an 
odd number of electrons. Reasonably high Kondo tem- 
peratures (i.e. Tk > 0.1 K) can be achieved when the QD 
is well coupled to normal- type leads. The Kondo effect 
manifests itself as a zero-bias conductance peak whose 
width is proportional to Tk- In 5'-QD-5' devices, the 
quasiparticle density of states (DOS) around the Fermi 



level (Ep) of the leads vanishes due to the opening of 
the superconducting gap (A). This lack of quasiparticles 
precludes Kondo screening. 

The competition between the Kondo effect and su- 
perconductivity is governed by the corresponding energy 
scales, UbTk and A. While no Kondo screening occurs 
for ksTK ^ A (weak coupling), a Kondo singlet is ex- 
pected to form for ksTK » A (strong coupling) at the 
expense of the breaking of Cooper pairs at the Fermi level 
[2^ . A quantum phase transition is predicted to take 
place at UbTk ~ A |13 - [l4| . Experimental signatures of 
this exotic crossover have been investigated both in the 
Josephson supercurrent regime [l8l - [20l and in the dissi- 
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pative sub-gap transport regime Yet a full un- 

derstanding of these experimental findings is still lacking. 
In this Letter, we report an experimental study on S-QD- 
S devices where the relative strength between Kondo and 
superconducting pairing correlations is tuned by means 
of a magnetic field, B, acting on the superconducting 
gap. The transition from strong to weak coupling is con- 
tinuously accomplished by sweeping B from just above 
the critical field. Be, where A = 0, to zero field, where A 
attains the maximum value, Aq, exceeding ksTK- 

The S-QD-S devices were fabricated from individual 
InAs/InP core/shell NWs grown by thermal evaporation 
(total diameter ^ 30 nm). The InP shell (thickness ^ 2 
nm) acts as a confinement barrier resulting in an en- 
hanced mobility of the one-dimensional electron gas in 
the InAs NW core After growth, the NWs were de- 
posited onto a degenerately doped, p-type Si substrate 
(to be used as a back gate), covered by a 300-nm-thick 
thermal oxide. Device fabrication was accomplished by e- 
beam lithography, Ar+ bombardment (to remove the NW 
surface oxide), metal evaporation, and lift-off. Source 
and drain contacts consisted of Ti (2.5 nm)/Al (45 nm) 
bilayers with a lateral separation of approximately 200 
nm, and a superconducting critical temperature of ?» 1 
K. Transport measurements were performed in a He3-He4 
dilution refrigerator with a base temperature of 15 mK. 
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FIG. 1. (color online) (a) Left panel: Color plot of dl /dV vs 
(B^, Vsd) measured at the center of diamond a. The dashed 
lines highlight the emergence of finite-bias peaks related to 
the opening of a superconducting gap. The superimposed 
line trace shows the B^-dependence of the linear conduc- 
tance. Right panel: dl / dV {Vsd) traces taken at different 
values. The inset shows a scanning electron micrograph 
of a typical device (scale bar: 200 nm). (b) dl /dViVsD) 
traces measured in an even diamond revealing the Dynes-like 
DOS of the strongly coupled lead, (c) Numerical calculations. 
The smaller solid dots denote the position of the A and 2A 
peaks, whereas the bigger open dots highlight the position of 
the Shiba bound state peaks, (d) Linear conductance (Go) 
normalized to the normal-state value {Gq ) and plotted as a 
function of A/ksTK for diamonds a (red dots) and /3 (black 
squares) . 



At low temperature, electron transport is dominated 
by the Coulomb blockade effect with the NW channel 
behaving as a single QD. Charge stability measurements 
(i.e. differential conductance, dl/dV, as a function of 
source-drain bias, Vsd, and back-gate voltage, Vc) were 
performed to identify Kondo resonances in Coulomb dia- 
monds with odd occupation and a spin- 1/2 ground state. 
These measurements were taken at 15 mK with the leads 
in the normal state (superconductivity was suppressed 
by means of a magnetic field B-^ = 70 mT perpendic- 
ular to the substrate and exceeding the perpendicular 



critical field B;^). In each Kondo diamond, Tk was mea- 
sured from the half width at half maximum (HWHM) 
of the zero-bias dl/dV peak, while the tunnel coupling 
asymmetry was extracted from the peak height, i.e. the 
linear conductance G, according to the relation: G/ 
Go = 'iTLTB/iTL + ri)2, where Go - 2eyh and Tl^r) 
is the tunnel coupling to the left (right) lead. Here we 
present data corresponding to three Kondo diamonds la- 
beled as a, (3 and 7, where: Tk a ~ 0.56 K, {TL/^R)a ~ 
6.6 X 10-3, Tk,p « 1 K, [Tl/Tr)p « 0.44, and Tk,^ « 
0.71 K, {Tl/Tr)^^1.&x 10-3. 

Figure la shows a dl /dV{B^ , Vsd) measurement taken 
at the center of Kondo diamond a. The zero-bias Kondo 
peak is apparent above B^ « 23 mT (we note that at 
such low fields the Zeeman splitting is much smaller than 
/cbTr-^q, explaining the absence of a split Kondo peak). 
Surprisingly, reducing the field below B^ does not lead to 
an abrupt suppression of the Kondo peak. Instead, the 
peak becomes progressively narrower and smaller, van- 
ishing completely only below B^ « 9 mT. 

We argue that the observed zero-bias peak is a manifes- 
tation of Kondo screening due to intra-gap quasiparticle 
states. To support this interpretation, we show in Fig. 
lb a data set taken in an adjacent diamond with even 
occupation (i.e. with no Kondo effect). The dI/dV{Vsd) 
traces shown correspond to different in- plane fields, S", 
ranging from zero to just above the in-plane critical field 

Be ~ 200 mT. When lowering the fields from above to 
below Be, the sub-gap dl/dV does not drop abruptly, 
supporting our hypothesis of a sizable quasiparticle DOS 
at Ep (we note that, although the measurement of Fig. 
lb refers to in-plane fields, a qualitatively similar behav- 
ior can be expected for perpendicular fields). The devel- 
opment of a "soft" gap just below the critical field is ad- 
ditionally marked by the absence of dl/dV peaks charac- 
teristic of the BCS DOS singularities. Such peaks develop 
only at fields well below _b| , becoming most pronounced 
at i? = 0. In this low-field limit, the sub-gap conduc- 
tance simultaneously vanishes and first-order multiple- 
Andreev-reflection resonances emerge at eVsd ~ ±A. 

To reproduce the observed sub-gap features and the co- 
existence of a Kondo peak and a superconducting gap, we 
calculated the dl/dV of a QD modeled by an Anderson 
Hamiltonian including coupling to BCS-type supercon- 
ducting reservoirs. We used the so-called non-crossing 
approximation, a fully non-perturbative theory that in- 
cludes both thermal and quantum fluctuations, comple- 
mented with the Keldysh-Green's function method to 
take into account non-equilibrium effects at finite Vsd (see 
Supplemental Material). In order to fit the experimen- 
tal data, the DOS of the leads was modeled by a Dynes 
function: 
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where "f{B) is a phenomenological broadening term. For 
small A(i?), the ^(B) term is particularly important 
leading to a finite quasiparticle DOS at Ep. By con- 
trast, as A{B) increases (with decreasing B), the DOS 
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of the superconducting leads approaches the ideal (i.e. 
unbroadened) BCS profile. 

Figure Ic shows a set of calculated dI/dV{Vsd) traces 
at different magnetic fields. By adjusting the DOS pa- 
rameters 7(i3) and A(i?), these calculations clearly show 
a zero-bias peak persisting below the critical field, in 
agreement with the experimental data of Fig. la. Since 
this peak emerges only in the case of a finite ^{B), we 
conclude that the finite DOS at the Fermi level is at the 
origin of the experimentally observed zero-bias anomaly. 
The narrowing of this Kondo anomaly with increasing A 
can be interpreted as a decreasing Tr- due to the shrink- 
ing quasiparticle DOS around the Fermi level (this as- 
pect is more quantitatively discussed further below). As 
Tk approaches the electronic temperature also the peak 
height gets smaller leading to the disappearance of the 
zero-bias peak. 

Figure 2 shows a second data set taken in Kondo di- 
amond /?. In this case, the stronger and more symmet- 
ric coupling to the leads results in a higher Tk, and a 
larger peak conductance. Nevertheless, the field depen- 
dence (Fig. 2a) shows substantially the same behavior 
as in Fig. la, i.e. a zero-bias Kondo peak persisting be- 
low , becoming progressively narrower with decreas- 
ing B-^, and vanishing below B^ w 10 mT. Interestingly, 
a sharp dl/dV resonance is found around Vsd = su- 
perimposed to the (wider) zero-bias Kondo peak. This 
resonance persists throughout the entire field range in 
which the leads are superconducting (see Figs. 2a and 
2b). By performing current-bias measurements (inset of 
Fig. 2b), we were able to ascribe this sharp resonance to 
a dissipationless Josephson supercurrent as high as 0.9 
nA at i? = 0. This finding reveals the possibility of a 
coexistence between the Josephson effect, linked to the 
superconducting nature of the leads, and a Kondo effect 
arising from the exchange coupling between the localized 
electron and intra-gap quasiparticle states. 

Temperature dependence measurements performed in 
diamond (3 (Figs. 2c and 2d, see also Supplemen- 
tal Material) show that in the regime of coexistence 
between the Kondo effect and superconductivity, the 
height of the persisting Kondo anomaly follows the 
usual T-dependence for a normal state regime given 
by G {T)/Go = [1/(1 + (T/T;,)']^ where = 
T*j^l\/2h - 1^ = 0.22 and is the effective Kondo 
temperature [2J] . In addition, we show that decreases 
with decreasing B^ (inset of Fig. 2d), as expected from 
the reduced intra-gap quasiparticle DOS. 

At base temperature, the peak heights of Kondo reso- 
nances a and /3 appear to follow approximately the same 
dependence on /^./ksTK (Fig. Id), where Tk refers to 
the Kondo temperature in the normal state. A similar 
scaling was reported earlier by Buizert et al. [l^, for 
a S-QD-S device fabricated from an InAs self-assembled 
quantum dot using Ti/Al contacts. In that paper, it 
was speculated that when fc^Tx 3> A, it becomes ener- 
getically favorable for Cooper pairs to split in order to 
screen the local spin and create a Kondo resonance at the 
Fermi level. The results presented here point at a differ- 



ent interpretation based on the presence of the already 
discussed intra-gap quasiparticle states, which become 
particularly important when B approaches Be- Accord- 
ing to this interpretation, the apparent scaling in Fig. 
Id is intimately related to a quasi-particle "poisoning" 
of the superconducting gap. 

Well below the critical field, as the zero-bias Kondo 
anomaly disappears, the dl / dV {Vsd) is dominated by a 
pair of peaks symmetrically positioned with respect to 
Vsd = (Fig. la). These peaks become most pronounced 
at i? = 0. Similar types of sub- gap structures have been 
reported in earlier works and they were given different 
interpretations: a Kondo enhancement of the first order 
multiple Andreev reflection process [isl . [l7| . or, in the 
case of asymmetrically coupled S-QD-S devices, a per- 
sisting Kondo resonance, created by the strongly coupled 
lead, which is probed by the BCS DOS of the second, 
weakly-coupled lead ^16']. Both of the above interpreta- 
tions invoke Kondo correlations to explain the observed 
sub-gap structure, even though, as we have pointed out, 
these correlations get suppressed as the superconducting 
gap reaches its largest value at i? = 0. More recently, 
it was shown that the finite-bias sub-gap structure can 
be explained in terms of tunneling through magnetic An- 
dreev bound states 25, 26]. 

The problem of a magnetic impurity interacting with 
a superconductor was addressed already in the sixties 
by Yu, Shiba and Rusinov [27] . It was shown that intra- 
gap Andreev bound states, now often called Shiba states, 
emerge as a result of the exchange coupling, J, between 
the impurity and the superconductor. As later confirmed 
by experiments based on scanning tunneling spectroscopy 
[28j . Shiba states emerge as pairs of peaks in the local 
DOS symmetrically positioned at energies -izEs relative 
to the Fermi level, where Eb depends on J and \Eb | < A. 

The zero-field dl / dV [Vsd) trace in the right panel of 
Fig. la exhibits a pair of dl/dV peaks at eVsd ~ ±1.4A, 
followed by negative dl/dV regions. Taking into account 
the strong asymmetry in the tunnel couplings, these fea- 
tures can be well explained in terms of a pair of Shiba 
levels with Eb = 0.4A, created by the strongly coupled 
S lead, and tunnel-probed by the weakly coupled S lead 
(see Fig. 3a). The observed dl/dV peaks result from 
the alignment of these Shiba levels with the BCS gap- 
edge singularities. Precisely, one dl/dV peak is due to 
the onset of electron tunneling from the Shiba level be- 
low Ep to the empty quasi-particle band of the S probe. 
The other peak is due to the onset of electron tunneling 
from the occupied quasi-particle band of the S probe to 
the Shiba level above Ep. Increasing \Vsd\ beyond these 
resonance conditions leads to a reduced tunneling prob- 
ability and hence a negative dl/dV. The Shiba-related 
features observed if Fig. la are very well reproduced by 
the numerical results in Fig. Ic. A significantly different 
sub-gap structure is found in the case of relatively low 
contact asymmetry, as clearly shown in the lower panel 
of Fig. 2b. Here, both S leads interact with the QD spin 
resulting in a stronger exchange coupling. We find a pair 
of dl/dV peaks at eVsd ~ ±A, which implies Eb ~ 0. 
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FIG. 2. (color online) (a) Out-of-plane field dependence mea- 
sured in diamond fi. (b) dl /dV line profiles taken at different 
magnetic fields. The inset is a voltage-current measurement 
carried out at B = 0, which reveal transport of a dissipation- 
less supercurrent in the device, (c) Evolution of the dl/dV 
measured at ~ 15 mT, as a function of the tempera- 
ture (T). The zero-bias Kondo peak and the gap-related side 
peaks show distinct T dependences. The former is strongly 
suppressed with increasing T. (d) Temperature dependence 
of the normalized conductance of the zero-bias Kondo peak 
measured either above (squares) or below (circles and trian- 
gles) Be'". G(0) denotes the peak height measured at base 
temperature. The inset reveals that the effective Kondo tem- 
perature Tff decreases with increasing A. 



In addition, we observe a rather complex set of smaller 
peaks most likely due to multiple Andreev reflection pro- 
cesses. 

In the weak-coupling limit, the energy of the Shiba 
states is related to the superconducting gap through 
Eb = A(l — x)/{l +x), where x = 3(7ri^_F and is 
the Fermi velocity [26]. Since J is prestmiably indepen- 
dent of B, the energy of the Shiba states should evolve 
proportionally to A(i?) as B is varied. We have verified 
this dependence with a measurement performed in dia- 
mond 7, where, as in diamond a, tunnel couplings are 
strongly asymmetric. Differently from diamond a, how- 
ever, the low-energy transport is characterized by two 
pairs of intra-gap dl/dV peaks, as shown by the zero- 
field curve in Fig. 3c. As in Fig. la, the most promi- 
nent peaks (at eVsd ^ ±1.4A) correspond to the align- 
ment of the Shiba levels with the BCS coherence peaks 
of the weakly coupled contact, and they are consistently 
followed by negative dl/dV dips. The weaker peaks at 
e-Vsd ~ ±0.55A (in fact only one of them is clearly vis- 
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FIG. 3. (color online) (a) Schematics of the formation of YSR 
bound states resulting from the interaction of the QD with the 
strongly-coupled lead, (b) Differential conductance measured 
in diamond 7 as a function of the in-plane magnetic field, (c) 
The solid line depicts the dl /dV taken at i? = 0. Two pairs 
of peaks are observed at |eVsd| = A -I- Eb,Eb- The dashed 
line shows the equivalent dl/dV taken in an even diamond, 
for comparison, (d) Numerical calculation of the dl /dV for B 
— 0. (e) Position of the intra-gap peaks as a function of the 
magnetic field. The inset shows the data before performing 
the rescaling in units of A. The dashed line displays the field- 
dependence of A. 



ible) can be interpreted as 'replicas' of the Shiba peaks. 
Such 'replicas' are expected when the Shiba levels line 
up with the Fermi level of the weakly coupled S lead (for 
eVsd — ^Eb), provided a non negligible density of quasi- 
particles is present throughout its superconducting gap. 
This is apparent from the calculated dI/dV{Vsd) trace in 
Fig. 3d, which agrees fairly well with the experimental 
one. 

The Shiba peaks and their replicas shift towards Vsd — 
as B II is increased. Their voltage positions are plotted 
in the inset of Fig. 3e as solid and open dots, respec- 
tively. For comparison, the A(i?ll) dependence measured 
in a non-Kondo diamond (with even electron number) 
is also plotted in the same inset (dashed line). Up to 
S'l 120 mT, the Shiba peaks and their replicas evolve 
proportionally to A(i?ll) in agreement with the theoret- 
ical prediction. Above i?" ^ 120 mT, this behavior be- 
gins to be affected by the increasing Zeeman splitting 
(w 0.3 eV/T) of the QD spin doublet. The main Shiba 
peaks get strongly suppressed and they seem to even- 
tually merge into the normal-state, Zeeman-split Kondo 
peaks through a non-trivial transition region (roughly be- 
tween 120 and 180 mT). In this parallel- field configura- 
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tion, the large Zeeman splitting at i?" w Be prevents the 
observation of a zero-bias Kondo peak when approaching 

In conclusion, we have studied the transport properties 
of a spin- 1/2 QD coupled to S contacts. The ability to 
continuously tune the superconducting gap with an ex- 
ternal magnetic field enabled us to investigate the tran- 
sition from a normal-state Kondo to a superconducting- 
state Shiba ground state. We showed that the presence 
of a finite quasiparticle DOS within the superconducting 
gap can promote the formation of a zero-bias Kondo peak 
coexisting with superconductivity. The origin of this 
quasiparticle DOS remains to be clarified, also through 
a better understanding of the superconducting proxim- 
ity effect in semiconductor nanowire structures (e.g., the 
role of disorder-induced pair breaking [1^). Finally, we 
should like to emphasize that our results bear clear impli- 



cations in the interpretation of sub-gap transport features 
in hybrid superconductor-semiconductor systems, espe- 
cially in the presence of relatively high magnetic fields 
that can cause a significant suppression of the supercon- 
ducting gap. This is precisely the regime where Majorana 
fermions are expected to arise as zero-energy quasiparti- 
cle states. We argue that intra-gap quasiparticle states, 
manifesting through a sizable background conductance 
(which is the case of experiments reported in Refs. 
could result in the screening of a local spin (or orbital) 
degeneracy leading to zero-bias anomalies that are not 
related to Majorana physics. 
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Science and Innovation through grant FIS2009-08744. 
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Model 



The dot is described by an Anderson Hamiltonian 



in ; Hd= e„dld„ + Un^ni , (1) 



where Sa is the single-particle energy level of the localized state with spin a, dj. (da-) the fermion creation (annihilation) 



. operator of the state, = d\dr, the occupation, and U the on-site Coulomb interaction, which defines the charging 

■ energy. We focus on a regime where the charging energy is much bigger than other energy scales. In this regime the 
^) . Hamiltonian in Eq. ([T|) suffices to describe all relevant physics. 

^ ■ Kondo physics arises as a result of the interplay between the strong correlation in the dot and the coupling of the 

^ ; . localized electrons with the electrons in conduction bands. In our case, these are described as two leads [a = L and 

' R) which can be either normal or superconducting. In this later case, they are represented by a BCS hamiltonian of 

^ . the form: 

o. 

^ O ^ with A being the superconducting pairing gap. Tunneling is described by the Hamiltonian 

^ : Ht = J2 {Vk^-<^Lada + h.c}j . (3) 

in . The total Hamiltonian is then given hy H = + Hn + Ht + . For simplicity, we ignore the k- and cr-dependence 

■ of the tunneling amplitudes. Therefore, we consider a simplified model with Vk^a — Va/\/2 which defines the widths 
' [ . = TrpolKip, where po is the (normal) density of states in the reservoirs. In the presence of superconductivity, the 

■ coupling to the reservoirs is modified due to the BCS density of states as 

O" 

^" V^{E)^T^N,{E,B)=T^Re ^ (4) 

E\'^- I\{Bf 

. This, of course, corresponds to a clean BCS case. As discussed in the main text, it is clear from the experimental data 

\^ ' in even diamonds that the clean BCS case is not a good description of the experiments. Instead, a Dynes expression 

. 5^ . of the form 

has to be used, where 7(i3) is a phenomenological broadening which takes into account a finite density of states inside 
the BCS gap. 



Non-crossing approximation method 

Now we write the physical fermionic operator as a combination of a pseudofermion and a boson operator as follows: 
d„ — b'^ fa where /.o- is the pseudofermion which annihilates one "occupied state" with spin a, and is a boson 



2 



operator which creates an " empty state" . We are interested in a hmit where the Coulomb interaction is very large 
such that we can safely take the limit of J7 — >■ oo. This fact enforces the constraint J^a fafa-+b''b = 1, that prevents the 
accommodation of two electrons at the same time in QD level. This constraint is treated with a Lagrange multiplier. 



+ A(^/t^ + 5t6-l) . (6) 



Notice that we have rescaled the tunneling amplitudes Va — > V aV N according to the spirit of a 1/A^-expansion {N 
is the total degeneracy of the localized orbital) . 

Our next task is to solve this Haniiltonian, which is rather complicated due to the presence of the three operators in 
the tunneling part and the constrain. Moreover, we need to take into account superconductivity and non-equilibrium 
effects. In order to do this we employ the so-called Non-Crossing approximation (NCA) fl'-'3| generalized to the 
superconducting case 0, Hj. Without entering into much detail of the theory, we just mention that the boson fields 
in Eq. ([6|) are treated as fluctuating operators such that both thermal and charge fluctuations are included in a self- 
consistent manner. In particular, one has to derive self-consistent equations-of-motion for the time-ordered double-time 
Green's function (sub-indexes are omitted here): 

iG{tX) ^ {TJ{t)p{t')), 

iB(t,t') ^ {TMt)b\t')), (7) 

or in terms of their analytic pieces: 

iG{t,t') = G>{t,t')e{t-t')-G<{t,t')0{t' -t), 

iB{t,t') = B>{t,t')d{t-t')+B<{t,t')d{t' -t); (8) 

A rigorous and well established way to derive these equations-of-motion was first introduced by Kadanoff and Bayni 
0, and has been related to other non-equilibrium methods (like the Keldysh method) by Langreth, see Ref. 8] for 
a review. In the paper, we just show numerical results of the coupled set of integral NCA equations for our problem 
and refer the interested reader to Refs. for details. In particular, the density of states is given by 

piu:)^--J2lm[A:iu;)], (9) 

TT ^ — ' 

(T 

where A'^^ie) is the Fourier transform of the retarded Green's function Al{t) = Gl{t)B<{~t) - G<{t)B''{-t). Note 
that this decoupling neglects vertex corrections and, as a result, the NCA fails in describing the low-energy Fermi- 
liquid regime. Nevertheless, the NCA has proven to give reliable results even at temperatures well below the Kondo 
temperature (of the order of T = IO^^Tk) 0- Following Meir and Wingreen in Ref. the current is given by: 

Io.^{L,R} = -^E/ deT^ie)[2ImAUe)Uie) + A<ie)]. (10) 
withA<(e) the Fourier transform of A<(t) = iG<{t)[B''{~t)-B°-{~t)] and/a(e) = the Fermi-Dirac function 

l+e kr" 

at each reservoir held at a chemical potential fia such that the applied bias voltage is defined as eV — — ■ 

In practice, we selfconsistently solve the NCA integral equations until good numerical convergence is reached. All 
dl/dV calculations presented in the main text are done for finite temperatures T = 0.25Tk and increasing values of 
A, roughly from A = (bottom curve of Fig. Ic) to A = 20Tk (top curve of Fig. Ic). 



"Soft" gaps in out-of-plane and in-plane magnetic fields 

Here we present measurements in which the magnetic field dependence of the superconducting gap A was probed 
both in out-of-plane and in-plane configurations. For this purpose, we have performed measurements using S-QD-N 
devices, which are identical to the devices discussed in the main text, except for the fact that one of the S aluminium 



3 



B-"{mT) B"(mT) 




Vsd (|jV) Vsd ([iV) 



FIG. 1: Magnetic field dependence of the superconducting gap measured in an out-of-plane (left panel) or in an in-plane (right 
panel) configuration. 
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FIG. 2: (a) Persistance and narrowing of the Kondo peak below Be. The side peaks observed in the red trace are related to 
the superconducting gap (see main text), (b) Temperature dependence of the persisting zero- bias Kondo anomaly, (c) dl /dV 
of the zero-bias peak plotted as a function of the temperature. 

leads is replaced by a normal metal Au contact. The data shown in Fig. 1 was taken in an even valley, i.e. in the 
absence of Kondo correlations. In agreement with our discussion in the main text, the transition from the normal 
to the superconducting regimes is characterized by the opening of a "soft" gap, which is populated by intra-gap 
quasiparticle states. As B decreases, the intra-gap DOS decreases, until a well-defined gap is obtained for low fields. 
The out-of-plane and in-plane field dependences are qualitatively identical, except for the different critical fields. 

Temperature dependence of a different Kondo resonance 

Here we present data corresponding to the temperature (T) dependence of a different Kondo resonance 5 {Tk.s ~ 
0.35 K). First, we verified that the overall magnetic field behavior of 5 is qualitatively identical to that observed in 
diamonds a, /3 and 7, as described in the main text. Fig. 2a depicts the persistance and narrowing of the zero- 
bias Kondo peak below the critical field Bjr . Side peaks related to the superconducting gap are also visible. The 
temperature dependence measurements were taken at = 10 mT. Fig. 2b clearly shows the distinct temperature 
dependences of the zero-bias peak and the side peaks. Indeed the zero-bias peak is strongly suppressed whereas the 
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side peaks are only weakly affected. This relates to the fact that the zero-bias peak follows a Kondo temperature 
dependence, which is expected to yield a 40% drop in conductance for T = 0.2 K and Tk,s ~ 0.35 K. The discrepancy in 
comparison with the experimental data (shown in Fig. 2c) most likely arises due to the the intra-gap quasiparticle DOS 
(see main text), which gives rise to a finite dl/dV at zero-bias even when the Kondo peak is thermally suppressed. 
On the other hand, the superconducting gap at i? = 0, Aq « 150 /ieV, is only reduced in « 10% at T = 0.2 K. 
Therefore, the position of the side peaks observed in Fig. 2b are only marginally affected by T. 
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